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Abstract 

A family of non-equilibrium statistical operators (NSO) is introduced 
which differ by the system lifetime distribution over which the quasi- 
equilibrium (relevant) distribution is averaged. This changes the form of 
the source in the Liouville equation, as well as the expressions for the 
kinetic coefficients, average fluxes, and kinetic equations obtained with 
use of NSO. It is possible to choose a class of lifetime distributions for 
which thermodynamic limiting transition and to tend to infinity of av- 
erage lifetime of system is reduced to the result received at exponential 
distribution for lifetime, used by Zubarev. However there is also other ex- 
tensive class of realistic distributions of lifetime of system for which and 
after to approach to infinity of average lifetime of system non-equilibrium 
properties essentially change. For some distributions the effect of "finite 
memory" when only the limited interval of the past influence on be- 
haviour of system is observed. It is shown, how it is possible to spend 
specification the description of effects of memory within the limits of 
NSO method, more detailed account of influence on evolution of system 
of quickly varying variables through the specified and expanded form of 
density of function of distribution of lifetime. The account of character of 
history of the system, features of its behaviour in the past, can have sub- 
stantial influence on non-equilibrium conduct of the system in a present 
moment time. 
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1 Introduction 



One of the most fruitful and successful ways of development of the description of 
the non-equilibrium phenomena are served by a method of the non-equilibrium 
statistical operator (NSO) [T] [3] [3J HI [5] . In work [5] new interpretation of a 
method of the NSO is given, in which operation of taking of invariant part 
PQ H HH [5] or use auxiliary "weight function" (in terminology 015]) in NSO 
are treated as averaging of quasi-equilibrium statistical operator on distribu- 
tion of past lifetime of system. In work [2] is chosen uniform distribution for an 
initial moment to, which after change of integration passes to the exponential 
distribution p q (u) — exp{— eu}. Such distribution serves as the limiting dis- 
tribution of lifetime 9J, of the first passage time of level. In general case it is 
possible to choose a lot of functions for the obvious type of distribution p q (u), 
that marked in works [6] [10]. This approach adjust with the operations spent 
in the general theory of random processes, in the renewal theory, and also with 
the lead Zubarev in work as [5] reception NSO by means of averaging on the 
initial moment of time. 

In Kirkwood's works [TT] it was noticed, that the system state in time 
present situation depends on all previous evolution of the non-equilibrium pro- 
cesses developing it. For example, in real crystals it is held in remembrance 
their formation in various sorts "defects", reflected in structure of the crystals. 
Changing conditions of formation of crystals, we can change their properties 
and create new materials. In works [5J it is specified, that it is possible to 
use many "weight functions" . Any form of density of lifetime distribution gives 
a chance to write down a source of general view in dynamic Liouville equation 
which thus becomes, specified Boltzmann and Prigogine 0[5J[I2], and contains 
dissipative items. 

If in Zubarev's works 0H] the linear form of a source corresponding limiting 
exponential distribution for lifetime is used other expressions for density of 
lifetime distribution give fuller and exact analogues of " integrals of collisions" . 

In works [13l [14] sources in the Liouville equation different from sources, 
entered in the NSO method in works [UEj are considered. But in works [13] [T4] 
at e — > (and r — > oo in interpretation paper [6]) this source turn into a zero. 
In the present work the sources which are not turn into a zero at T — > oo 
are considered. Such source is caused by external influences on system. The 
approach of present paper differs from the methods used in [13] [14] . But use of 
distribution of a lifetime of system in the present work can be compared with 
noted in work [T3] enlarged the set of macroobservables, include besides the 
physically necessary macroobservables additional ones, namely lifetime. 

In work [15] irreversible transfer equations are received in assumptions of 
coarsening of the distributions, a certain choice of macroscopical variables and 
the analysis of division of time scales of the description (last circumstance 
was marked in |16j). Importance and necessity of the analysis of the time 
scales playing a fundamental role in the description of macroscopical dynamics 
of system is underlined. Evolution of slow degrees of freedom is described 
by Markovian equations. Thus the time scale on which observable variables 
evolve, should be much more time of memory on which the residual effects 
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brought by irrelevant degrees of freedom are considered. As marked in paper 
[14] , it requires the condition of a maximum damping of the non- macroscopic 
information. 

Otherwise effects of memory play an essential role. Memory time is esti- 
mated in work [15] for Boltzmann equation. In the present work the consider- 
ation subject is made by situations when it is necessary to consider effects of 
memory. Examples of such situations are given in [15] . 

In work [10] it is shown, in what consequences for non-equilibrium properties 
of system results change of lifetime distribution of system for systems of the 
limited volume with finite lifetime. In the present work are considered also 
infinitely greater systems with infinite average lifetime. 

2 New interpretation of NSO 

In [5] the Nonequilibrium Statistical Operator introduced by Zubarev [TJ [5] 
rewritten as 



/•OO 

lng(t) = / p q (u)lng q (t - u, ~u)du, lng q (t,0) = - V"F n (t)P„; 

ln Qq {tM) = e^ H / m hng q (t,0)e^ H /^; *(t) = InSpexp^F^P^, 

n 

(i) 

where H is hamiltonian, lng(t) is the logarithm of the NSO in Zubarev's form, 
lng q (t,0) is the logarithm of the quasi-equilibrium (or relevant); the first time 
argument indicates the time dependence of the values of the thermodynamic 
parameters F m ; the second time argument ti in g q (t\,t2) denotes the time 
dependence through the Heizenberg representation for dynamical variables P m 
from which g q (t, 0) can depend [TJ[21[B]|71[5]. In [B] the auxiliary weight function 
p q (u) = eexp{— eu) was interpreted as the density of probability distribution 
of lifetime T of a system. There Y is random variables of lifetime from the 
moment to of its birth till the current moment t\ e^ 1 = (t — to); (t — to) = (Y), 
where (r) = J up q {u)du is average lifetime of the system. This time period 
can be called the time period of getting information about system from its 
past. Instead of the exponential distribution p q (u) in (1) any other sample 
distribution could be taken. This fact was marked in [5] and [71 [5] (where the 
distribution density p q (u) is called auxiliary weight function w(t, t')). From the 
complete group of solutions of Liouville equation (symmetric in time) the subset 
of retarded "unilateral" in time solutions is selected by means of introducing a 
source in the Liouville equation for In g(t) 

+ ih In g{t) = - £ (ln g(t) - In g q {t, 0)) - J, (2) 

which tends to zero (value e — > 0) after thermodynamic limiting transition. 
Here L is Liouville operator; iL = -{H,g} = E fc [|^^ - f§^-]; H is 
Hamilton function, p^ and are momenta and coordinates of particles; {...} is 
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Poisson bracket. In [T7] it was noted that the role of the form of the source term 
in the Liouville equation in NSO method has never been investigated. In [18] it 
is stated that the exponential distribution is the only one which possesses the 
Markovian property of the absence of contagion, that is whatever is the actual 
age of a system, the remaining time does not depend on the past and has the 
same distribution as the lifetime itself. It is known [TJ [3J [5] that the Liou- 
ville equation for NSO contains the source J = J zu b — —s[lng(t) — lng q (t, 0)] 
which becomes vanishingly small after taking the thermodynamic limit and 
setting e — > 0, which in the spirit of the paper PQ corresponds to the infinitely 
large lifetime value of an infinitely large system. For a system with finite size 
this source is not equal to zero. In [5] this term enters the modified Liou- 
ville operator and coincides with the form of Liouville equation suggested by 
Prigogine |12] (the Boltzmann-Prigogine symmetry), when the irreversibility 
is entered in the theory on the microscopic level. We note that the form of 
NSO by Zubarev cast in [6] corresponds to the main idea of [12] in which one 
sets to the distribution function g {g q in Zubarev's approach) which evolves ac- 
cording to the classical mechanics laws, the coarse distribution function g (g(t) 
in the case of Zubarev's NSO) whose evolution is described probabilistically 
since one perform an averaging with the probability density p q (u). The same 
approach (but instead of the time averaging the spatial averaging was taken) 
was performed in [T5] , 

Besides the Zubarev's form of NSO [TJ[5], NSO Green-Mori form [HI [21] 
is known, where one assumes the auxiliary weight function [7] to be equal 
W(t,f) = 1 - (t - t')/T;w(t,t') = dW(t,t')/dt' = l/r;r = t - t . After 
averaging one sets r — > oo. This situation at p q {u = t — to) = w{t,t' = to) 
coincides with the uniform lifetime distribution. The source in the Liouville 
equation takes the form J = lng q /r. In [TJ this form of NSO is compared to 
the Zubarev's form. 

One could name many examples of explicit defining of the function p q (u). 
Every definition implies some specific form of the source term J in the Liouville 
equation, some specific form of the modified Liouville operator and NSO. Thus 
the family of NSO is defined. 

3 Modifications to the nonequilibrium descrip- 
tion 

Let's consider now, what consequences follow from such interpretation of NSO. 
3.1 Families of NSO 

Setting various distributions for past lifetime of the system, we receive a way of 
recording of families of NSO. Class of NSO from this family will be connected 
with a class of distributions for lifetime (taken, for example, from the stochastic 
theory of storage processes, the theory of queues etc.) and with relaxation 
properties of that class of physical systems which is investigated. The general 
expression for NSO with any distribution 
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lng(t) = / p q (u)lng q (t — u, —u)du = 
Jo 



(3) 



f°° i f . (iln — u, — u) 
= lng q {t, 0) - J (J p q {u)du) — du, 

where integration by parts in time is carried out at J p q (y)dy\ y=0 = — 1; 
J Pg(y)dy\ y ^ 00 = 0; at p q (y) = eexp{-ey};e = 1/(T), the expression (1) 
passes in NSO from [HE]- I n [H] it is shown, how from random process 
X(t), corresponding to evolution of quasi-equilibrium system, it is possible to 
construct set of new processes, introducing the randomized operational time. 
It is supposed, that to each value t > there corresponds a random value 
r(t) with the distribution p q (y). The new stochastic kernel of distribution 
of a random variable X(T(t)) is defined by equality of a kind (1). Random 
variables X(T(t)) form new random process which, generally speaking, need not 
to be of Markovian type any more. Each moment of time t of "frozen" quasi- 
equilibrium system is considered as a random variable T(t) the termination of 
lifetime with distribution p q (y). Any moment of lifetime can be with certain 
probability the last. That the interval t — to = y was enough large (that 
became insignificant details of an initial condition as dependence on the initial 
moment to is nonphysical [HE]), it is possible to introduce the minimal lifetime 
Tmm = Ti and to integrate in ([3]) on an interval (ri,oo). It results to the 
change of the normalization density of distribution p q (y). For example, the 
function p q {y) = eexp{— ey} will be replaced by p q (y) — eexpjeri — ey}, y > 
^i>Pq(y) — 0, y < Ti- The under limit of integration in (JSJ by Ti — > is equal 0. 
It is possible to choose p q (y) = Cf(y),y < h;p q {y) = e cxp{-ey}, y >h;C = 
(1 — exp{— eii})/(/ ' 1 f(y)dy). The function f(y) can be taken from models 
of the theory of queues, the stochastic theory of storage and other sources 
estimating the lifetime distribution for small times (for example [22[ [23| [9l [24] ) . 
The value t\ can be found from results of work [9]. It is possible to specify many 
concrete expressions for lifetime distribution of system, each of which possesses 
own advantages. To each of these expressions there corresponds own form of 
a source in Liouville equation for the noncquilibrium statistical operator. In 
general case any functions p q {u) the source is: 

J = p,(Q) lng q (t, 0) + £ ®EM(i nQq (t - y, - y ))dy (4) 

(when values p q (0) disperse, it is necessary to choose the under limit of inte- 
gration equal not to zero, and r m , n ). Such approach corresponds to the form 
of dynamic Liouville equation in the form of Boltzmann-Bogoliubov-Prigogine 
[7] [5] [H] , containing dissipative items. 

Thus the operations of taking of invariant part pQ, averaging on initial 
conditions [2j , temporary coarse-graining [TT] , choose of the direction of time 
[8], are replaced by averaging on lifetime distribution. 

It is essentially that and in exponential distribution from [HE] e ^ 0. The 
thermodynamic limiting transition is not performed, and actually important 
for many physical phenomena dependence on the size of system are considered. 
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We assume e and (r) to be finite values. Thus the Liouville equation for g(t) 
contains a finite source. The assumption about finitencss of lifetime breaks 
temporary symmetry. And such approach (introduction p q (y), averaging on it) 
can be considered as completing the description of works [TJ [2] ■ 

In work [5] lifetimes of system are considered as the achievement moments 
by the random process characterizing system, certain border, for example, zero. 
In [9] are received approached exponential expressions for density of probability 
of lifetime, accuracy of these expressions is estimated. In works [2"5l |2"o] life- 
times of molecules are investigated, the affinity of real distribution for lifetime 
and approached exponential model is shown. It is possible to specify and other 
works (for example [23 HH1 US]) where physical appendices of concept of life- 
time widely applied in such mathematical disciplines, as reliability theory, the 
theory of queues and so forth (under names non-failure operation time, the em- 
ployment period, etc.) are considered. In the present section lifetime joins in a 
circle of the general physical values, acting in an estimation or management role 
(in terminology of the theory of the information |30| ) for the quasi-equilibrium 
statistical operator that allows to receive the additional information on system. 

Let's notice, that in a case when value d\ng q (t — y, —y)/dy (the operator 
of entropy production a pQ) in the second item of the right part © does not 
depend from y and is taken out from under integral on y, this second item 
becomes a(T), and expression (|3J) does not depend on form of function p q {y). 
There is it, for example, at g q (t) ~ exp{— at}, a = const. In work [21] such 
distribution is received from a principle of a maximum of entropy at the set of 
average values of fluxes. 

3.2 Physical sense of distributions for past lifetime of sys- 
tem 

As is known (for example, [22], [9], [21]), exponential distribution for lifetime 

Pq(y) = eexp{-e?/}, (5) 

used in Zubarev's works [TJ [5], is limiting distribution for lifetime, fair for 
large times. It is marked in works [U [2] where necessity of use of large times 
connected with damping of nonphysical initial correlations. Thus, in works 
[U [2] the thermodynamic result limiting and universal is received, fair for all 
systems. It is true in a thermodynamic limit, for infinitely large systems. 
However real systems have the finite sizes. Therefore essential there is use of 
other, more exact distributions for lifetime. In this case the unambiguity of the 
description peculiar to a thermodynamic limit is lost |32j . 

For NSO with Zubarev's function (|5j the value enter in second item 

-/ P,(y)dy = exp{- £ y} = l-ey+(ey) 2 /2-... = 1-y / (T)+y 2 /2(F) 2 - ... (6) 

Obviously that to tend to infinity of average lifetime, (r) — » oo, correlation © 
tends to unity. 
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Besides exponential density of probability ([5]), as density of lifetime dis- 
tribution Erlang distributions (special or the general), gamma distributions 
etc. (see [22J 123]), and also the modifications considering subsequent com- 
posed asymptotic of the decomposition [53] can be used. General Erlang dis- 
tributions for n classes of ergotic states are fair for cases of phase transitions 
or bifurcations. For n — 2 general Erlang distribution looks like p q (y) — 
6pi exp{— piy} + (1 — 0)p2 exp{— piy}, < 1- Gamma distributions describe 
the systems which evolution has some stages (number of these stages coincides 
with gamma distribution order). Considering real- life stages in non-equilibrium 
systems (chaotic, kinetic, hydrodynamic, diffusive and so forth 16J ) , it is easy 
to agree, first, with necessity of use of gamma distributions of a kind 

Pq (y) = eieyf- 1 exp{- £ y}/r(fc) (7) 

(T(k) is gamma function, at k = 1 we receive distribution J^J), and, secondly, 
- with their importance in the description of non-equilibrium properties. The 
piecewise-continuous distributions corresponding to the different stages of evo- 
lution of the system will be used below. 

More detailed description p q (u) in comparison with limiting exponential 
(J5| allows to describe more in detail real stages of evolution of system (and 
also systems with small lifetimes) . Each from lifetime distributions has certain 
physical sense. In the theory of queues, for example [33 , to various disciplines of 
service there correspond various expressions for density of lifetime distribution. 
In the stochastic theory of storage [33], to these expressions there correspond 
various models of an exit and an input in system. 

The value e without taking into account of dissipative effects can be defined, 
for example, from results of work [9 . The value e is defined also in work [6] 
through average values of operators of entropy and entropy production, flows 
of entropy and their combination. 

How was already marked, it is possible to specify very much, no less than 
1000, expressions for the distributions of past lifetime of the system. Certain 
physical sense is given to each of these distributions. To some class functions of 
distributions, apparently, some class of the physical systems corresponds, the 
laws of relaxation in which answer this class of functions of distributions for 
lifetime. 

3.3 Influence of the past on non-equilibrium properties 

In |18j by consideration of the paradox connected with a waiting time, the 
following result is received: let X\ — Sx^Xi = S2 — Si;... are mutually in- 
dependent also it is equally exponential the distributed values with average 
1/e. Let t > is settled, but it is any. Element Xk, satisfying to condi- 
tion Sfc_i < t < Sk, has density v t (x) — e 2 xexp{— ex}, < x < t;v t (x) = 
e(l + ex) exp{— ex}, x > t. In Zubarev's NSO [TJ [2] the value of lifetime to a 
present moment t, belonging lifetime Xk, influence of the past on the present 
is considered. Therefore the value p q {u) should be chosen not in the form of 
exponential distribution (J5j, and in a form 
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Pq(y) = £ 2 yexp{~ey}, 



(8) 



that in the form of gamma distribution at k = 2. In this case distribution 
((5J) coincides with special Erlang distribution of order 2 [55], when refusal (in 
this case - the moment t) comes in the end of the second stage [22], the sys- 
tem past consists of two independent stages. Function of distribution is equal 
P q {x) = 1 — exp{— ex} — exexp{— ex},p q (u) = dP q (u) / du, unlike exponential 
distribution, when P q {x) = 1 — exp{— ex}. The behaviour of these two densities 
of distribution of a form ((SJ) and © essentially differs in a zero vicinity. In case 
of © at system low probability to be lost at small values y, unlike exponential 
distribution §5§ where this probability is maximal. Any system if has arisen, 
exists any minimal time, and it is reflected in distribution ((8]). 

In work 1 10] it is shown, in what differences from Zubarev's distribution ([3]) 
with exponential distribution of lifetime ([5]) results gamma distribution ([7]), 
([8j use. Additional items in NSO, in integral of collisions of the generalized 
kinetic equation, in expressions for average fluxes and self-diffusions coefficient 
are considered. The same in [10] is done and for special Erlang distribution k = 
2,3,4...,n,P q (x) = l-cxp{-ea;}[l + ex/l! + ... + (ea;) fc - 1 /(fc-l)!];e = k/(T). 
For distributions (O, © is correct correlation (4), value — J p q (u)du — » 1 by 



Thus the multi-stage model of the past of system is introduced. Non- 
equilibrium processes usually proceed in some stages, each of which is char- 
acterized by the time scale. In distribution ([8]) the account of two stages, 
possibly, their minimal possible number is made. Other distributions can de- 
scribe any other features of the past. Corresponding additives will be included 
into expressions for fluxes, integral of collisions, kinetic coefficients. Besides 
special Erlang distributions with whole and specified values k = n, that does 
not deduce us from set of one-parametrical distributions, the general already 
two-parametrical gamma distribution where the parameter k can accept any 
values can be used. In this case (r) = k/e. The situation (formally), when 
k < 1 is possible. Then sources will tends to infinity, as (t — fcOkZ^ ^ 00 
at k < 1. This divergence can be eliminated, having limited from below the 
value t — to of minimal lifetime value T m i n , having replaced the under zero limit 
of integration on T m i n - Then to expression for a source ((4]) it is added item 
[(eT min ) k - 1 /T(k)]e exp{-er mm } In g q (t - 

In [18] it was shown that the exponential lifetime distribution tf — to (tf , 
to are random moments of system death and birth) at big t's the "age" of a 
system t — to tends to the exponential form. In Zubarev's NSO [TJ [2] the lifetime 
value t — to to the current time t, which is a part of the total lifetime tf — t , 
is considered, that is the influence of the past on the current moment is taken 
into account. The full lifetime distribution, as well as the "past" lifetime (i.e. 
time from the system birth to till the current time t) need not be exponential. 
The logarithm of NSO in the case ([8]) has the form 



<r> -> oo. 




(9) 
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poo 

/ £ 2 uexp{— eu} In g q (t — u, —u)du = 
Jo 

/•OO 

lng q (t, 0) + / a(t — u, — u)(l + eu) exp{— eu}du ■ 



n 



\nQzub{t)+ a(t—u, —u)euexp{—eu}du; a(t—u,—u) = dlng q (t—u,—u)/du, 



where lng zu b(t) — lng q (t,0) + 

J °° cr(t—u, —u)exp{—eu}du = In g q {t—u 1 —u)e exp{— eu}du is the Zubarev's 
form of the NSO, a{t) = 9 In g q (t — u, —u)/du\ u= o — —d\nq(t,0)/dt is the en- 
tropy production operator pQ. It is seen from (jUJ) that the logarithm of the 
NSO has an additional term in comparison to the Zubarev's form. The source 
in the rhs of the Liouville equation (or dissipative part of the Liouville operator 
pfj) equal J = — e[lng(t) — lng zu b{t)}, that is the system relaxes not towards 
lng q (t,0), like it is the case of Zubarev's NSO, but towards lng zu b(t). From 
the expression © it is seen that introduced NSO contains amendments to the 
Zubarev's NSO [BUI- The physical results obtained with use of ([9]) also con- 
tains additional terms in comparison to Zubarev's NSO. The additional terms 
describe the influence of the lifetime finiteness on the kinetic processes. The 
expressions for average fluxes pQ averaged over © have the form 

<j m (x) >=<j m (x) > zub + (10) 

]T / T e(t~t')exp{e(t' -t)}(j m (x),j n (x',t' -t))X m (x',V)dt'dx', 

where < j m (x) > zub =< j m (x) > t + 

E„/ /_oo ex P{ £ (*' ~ t )}{j m { x )^3 n { x '^' ~ t)X m (x',t')dt'dx' are fluxes in the 
form obtained by Zubarev pQ, j n are flux operators, X m are corresponding 
thermodynamical forces; (j rn (x),j n (x',t)) — P^ 1 Jq < j m (x)(j n (x' ,t,ir) — 
< j n (x',t) >i) >i dr are quantum time correlation functions, j n (x' ,t,ir) — 
exp{-0- l AT}j n (x f ,t)exp{0- 1 AT};A = £ m / F m {x,t)P m {x)dx. The colli- 
sion integrals of the generalized kinetic equation pQ , averaged over © have the 
amendments 



S [ aL = -h- 2 / dtetexp{et} < [H x (t), [H u P k ] + ih £ P * ] > q 

J — oo m ° *m > 

(11) 

to Zubarev result pQ: S^ 2 ) = -/i~ 2 /"^ dtexp{et} < [H^t), [#i,P fc ] + 

dS (1> 

*^Em ^ m d<p > ] > «' wnere tne Hamiltonian of the system is H = H +Hi, Hi 
is the Hamiltonian of the interaction which contains the longtime correlations 
[7], = < ^ Pk ^ 1 ^ >q ■ The same is valid for the generalized transport equations 
PQ, kinetic coefficients etc. Thus the selfdiffusion coefficient (or, to be exact, 
its Laplace transform over time and space) obtained in [5] in the form 
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D(u,q) 



2 $(w,q) 



(12) 



[1 + 



(iu — e) 



where *( w , g ) = J "°° ^^"^ r^t'A"^^ » »■ = 



At e — > 0, for infinitely large system in the thermodynamic limit this ex- 
pression (fT3|) coincides with (fT2]) at e — > 2 . For finite size systems (as well 
as for the case iv ~* 0) the results differ. 



Above, as well as in work [TO], additives to NSO in the Zubarev's form are 

received for systems of the finite size, with finite lifetime. We will show, as for 
systems with infinite lifetime, for example, for systems of infinite volume, after 
thermodynamic limiting transition, the same effects, which essence in influence 
of the past of system, its histories, on its present non-equilibrium state are fair. 

In work |10j it is shown, as changes in function p q (u) influences on non- 
equilibrium descriptions of the system. But for those distributions p q (u) , which 
are considered in [TD] (gamma-distributions, (J7J), ©) the changes show up only 
for the systems of finite size with finite lifetimes. Additions to unit in equation 
([SJ becomes vanishingly small to tend to infinity of sizes of the system and its 
average lifetime, as in the model distribution ([5]) used in Zubarev's NSO ©• 
For the systems of finite size and the exponential distribution results to nonzero 
additions in expression ([3]). Thus, these additions to NSO and proper additions 
to kinetic equations, kinetic coefficients and other non-equilibrium descriptions 
of the system, are an effect finitcness of sizes and lifetime of the system, not 
choice of distribution of lifetime of the system. We will find out, whether there 
are distributions of lifetime of system for which and for the infinitely large 
systems with infinitely large lifetime an additional contribution to NSO differs 
from Zubarev'a NSO. 

We will consider a few examples of task of function p q (u) in ((2|), ([3]). We 
will be limited to the examples of task of the piecewise-continuous distributions, 
which result in results different from works [TJ [5] . There are numerous exper- 
imental confirmations of such change of distribution of lifetime of the system 
p q (u) on the time domain of life of the system. This and transition to chaos 
and transition from the laminar mode to turbulent one is also accompanied 
by the change of distribution of p q (u). In works [34|, 135] is shown transition of 
distribution of the first passage processes from Gaussian regime to the not Levy 
conduct in some point of time. Besides the real systems possess finite sizes and 
finite lifetime. Therefore influence of surroundings on them is always present. 




(13) 



4 Systems with infinite lifetime 
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That a source not is equal to the zero and in a limit infinitely large systems, is 
related to the openness of the system, influence on her of surroundings 
4a). We will put 

p q (u) = { p^f' U< \ (14) 
oeexpj— eu\, u>c; 

First part of distribution (fT4| . Pareto distribution, in work [23] got from 
the exponential distribution with the random parameter of intensity p, when 

/>oo 

M«) = / pe- pu f p (p)dp 
Jo 

at fp(p) = a ^ ka ^ r(k)^~ ap ^ 1S g amma -distribution, a = -j^, p is average of 
gamma-distribution. From the condition of the normalization J °° p q (u) = 1 
we will find a normalization multiplier b = e £C {-^^) k . Average lifetime (r) for 
distribution Q14[) is equal 

The value ([15]) (r) -> oo at e -> 0. From |14]) we find 

6^r-«.r«>c (i6) 

Source in right part of Liouville equation @ for distribution (TT4"]) in accor- 
dance with expression ^ equal 



J = - In g q (t, 0) - f ^ k + 1 2t sdu - e C£ ( — ) k [°° e 2 e~ £U Sdu, (17) 

where S = In g q (t — u, —u). Distribution of NSO ([3]) in the case of distribution 
equal 

, Z" 00 „. dhx p(t — u, —u) 
lng(t) = lng q (t,0) + e~ eu ^- - -du + A; (18) 

Jo 

A= / [(— — ) fc -e- en ]adw+ / [e ec (^—) k - l]e- £U adu, 
Jo a + u J c a + c 

where a = dyn e(t-u,-u) , i Uzub g^ — lng q (t, 0) + Jq 00 e~ eu adu is distribution got 
by Zubarev in [TJ[2], an( i A are amendments to him. 
4b). We will consider distribution of kind now 

P ^ = i teM-eu}, u>c; (19) 

From the condition of the normalization we find a value d = -(1 — e~ £C ). 
Average lifetime 

(T> = ^- + V" (1+ec). (20) 
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Average lifetime (r) — > oo at e — > 0. Source Q in equation (J2J) in this case 
equal 

/>oo 

J = d\ng q (t,0)- J e 2 e- eu Sdu. (21) 

NSO is equal 

lng(t) = lng zub (t,0)- [ [er 611 + -(1 - e - £C )u}adu; (22) 



[ e - ev + -(1 - e- ec )u]adu. 
c 



We get in this case, that in an additional item memory of the system is 
limited by a size, is observed effect of limited memory. It is possible to consider 
and other examples of task distributions p q {u) which reduce to limited memory. 

4c). For the exponential density of distribution with different intensities in 
different temporal intervals 

p q (u) = { «< c - (23) 

from the condition of the normalization finds that b = e c ( £2_El ); 

<r) = — [1 - e- £lC (l + e lC )] + -e- Elc (l + e 2 c). (24) 

El £2 

At (r) -> oo, e 2 -> 0. 

J = stlng g (t,0) -J e\e- £lU Sdu-e c(e2 - £l) J e^e _£aU 

Jnfi(t) - lng zub (t) + A; (25) 

/•C /"OO 

A = y [e- eiu - e- £2U ]adu + J [ e < e *- e J - l]e- E2U <rdu; 



\e- £lU - l\adu. 



4d). We will consider yet distribution of kind 

e 2 uexp{— eu}, u < c, 
P«(«) = { b cxp{ ~ 7 " } ,— , u>c. ( 26 ) 



The second part of distribution ([26)1 can be received from results of works [36] 

T°[ a 7"] » [l-exp{-ec}(l+exp{-£c»] 



and like to Tsallis distributions 1 37 . For ([261) 6 = = = j- T — rv- 

where B(l — p, 1) is incomplete beta function [38] . p = 07 exp{— 700}. 



Average lifetime 



(r> = -[1 - e- £C (l + ec + (ec) 2 /2)] + 
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(1 ' 



where T(.) is gamma function, n F m is hypergeometrical functions [38] , tends 
to infinity at 7 — > 0. In this case value of the normalization b — > 0. Therefore 
at (r) — > 00, 7 — > and 6^0; 



Zng(t) = Zng 9 (t, 0) + / e~ eu (l + eu)adu; J=- e~ eu e 2 {\ - eu)Sdu. 



The effect of limited memory shows up in this case obviously. 

Why the examples of this section differ from examples of section 3? In 
interpretation [5] fluctuate the random value to in u = t — to- In [2] lim- 
iting transition is conducted for the parameter e, e — * in the exponential 
distribution p q (u) = eexp{— eu} after thermodynamical limiting transition. In 
interpretation of work [6] it corresponds to that mean lifetime of the system 
(r) = (t — to) = 1/e — > 00. But middle intervals between random shoves in- 
finitely increase, exceeding lifetime of the system. Therefore an item with a 
source in Liouville equation applies in a zero. If there is the change of distri- 
bution p q (u) on the time domain of life of the system, as in examples 4a)-4d), 
influence of surroundings, followed this change with, remains in the time do- 
main of life even at tendency to infinity of mean lifetime. 

5 Application of distributions of section 4 to 
research of conductivity 

On the example of determination of conductivity we will consider, in what 
consequences the change of type of functions p q (u) and g(t) results as compared 
to an exponential law for p q (u), used in work [5]. 

Determination of conductivity by a method NSO is considered in works 
[39l l40l |4T1 [42] . In this section we will consider the transport of charges in 
the electric field, as linear reaction on mechanical perturbation, conductivity 
in the linear approaching, following results of work and, as in [5], we will be 
limited to the important special case - reaction of the equilibrium system on 
the spatially homogeneous variable field. 






Hamiltonian of perturbation is given by a formula 




P = ejfj, 
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where is charge of particle, and fi is its radius vector. The operator of 
current is equal 

j = p = e Y,^ = ^X^> 

3 3 

where pj is particle momentum, m is mass of a particle. We will choose a model 
in which coulomb interaction is taken into account as self-consistent screening 
the field, i.e. we will consider that E = E°. Most essential difference from 
work [5] at consideration of this problem consists of replacement of Laplace 
transformation used in [5], form 

/>oo 

(A;B) u+is = dte^+^ t (A(t),B),(e>0), (27) 



where (A(t),B(t')) = J* dxTr{AA(t)AB(t' +i/3hx)g eq } is temporal correlation 
function [5], by other integral transformation. So, for an example 4a) with 
distribution of form (fI3|). iflB]) expression (|2"7)l it is replaced on 

(A; B) w;a , fe + e £C (^-) fe (A; B) u+ie .( Ci0o) , (28) 
a + c 

where 

PC />oo 

(^;BU= / ^e- t (-^-) fc (^),S);(A; J B) w+fe . (C)00) = / dte^ +i ^ (A(t) . 
Jo a + t Jc 

(29) 

We will consider a isotropic environment in which tensor of conductivity is 
diagonal. In work [5] for Laplace transformation of kind (|27j) expression is got 
for specific resistance p(u>) of kind 

p{u) = -L = [-iw + M] ; (30) 

o-(w) (3(J,J) 

M= ( (31) 

(J, J) + (J; J) u+le 

where V is the volume of the system, (3 is inverse temperature, <j{oS) is scalar 
coefficient of conductivity or simply conductivity. In the examples considered 
below in expression (j3Tj) a value M will change. Doing the operations conducted 
in work [5] , with replacement of expression ([27)) by expression (|28| , we will get 
for this case in place of correlation (|3ip more difficult expression of kind 

(*^! J)u>;a,k + e '" C ( ) h (J', J) a)+ie;(c,oo) 

- ^ ; 

X = (J(0), J) - -{J ; j) w;o fe+i - (1 - -^J— )e- c (-^) fe (J(c), J)+ 
a i(w + ze) a + c 

+(■/; J/ui;a,k H v ; r^r6 \J] J) u+ie-.(c,co)- 
1{U1 + IE) 
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At e -> 

(J; f) US a,k + e iuc (^ c ) k (j; J)^(coo) 



M = 



(J(0), J) — -(J; J)u;a,k+1 + (J', J)u;a,k + (J', J)u;(c,oc) 

For distribution (TIT)]) in 4b) Laplace transformation of kind (p?T|) it is replaced 



on 



-d(A; B)u; C .t + (A; B) u+i£ .( CtOQ ), 



where (A; B) u . c j = dte tu >H(A(t) , B) , the value d is given in PD|) . 
(vl; -B) w+fe .( C)tx) ) is given in (J25|. In place of expression (pTTj) in this case we 
will get expression 

M i(J+ie) > ^)w+ie;(c : oo) ~ d[w ! «i)w;c,i + (J 5 ^)w;c,t=l] 



if] 



IU1 '-> — ^ -< 

i^l = "7 . ■ \{J'i J}u+ic;(c,oo) — d[(J] J}cu- C ,t + 

t(Lu + ie) 

+ (J; J)u;c,t=i] + (dc + — ^—e- ec )e iuc (J(c); J). 
i{u> + is) 



If e ->■ 0, (r) -> oo, d -> 0, and 

(J; J) 

oj;(c, oo) 



M 



e^ c (J(c),J) + (J; J) w;(c ,oo) 



This expression at the small values c low differs from (|3Tj) . For a case 4c) with 
distribution p q (it) of kind |23|) Laplace transformation (|27|) it is substituted by 
a value 

(^4;-B}w+ie;(o : c) + e c ( E2_ei )(^4; i?) w +i e2 ;(c,oo) ; 

where (A;B) u+ie ,^) = dte i ^ +i ^ t (A{t), B), (4;-B) u+fea .( Cl0O ) is given in 
(|29p and value M from (13 lj) it is substituted by a value 



if 2 

X 2 = (J(0), J) - e ^+ iE i> c (J( C ); J) + ( j; j) w+feli(0)C ) + 

+e c(e 3 - £l )^±i£ll [e i(-+^)c ( j (c) j) + ( j. J) u+fe . ( j]. 
t(w + Z£ 2 ) 

At (r) — > oo and £2 — ► the value M changes unessential, assuming a form 

_ (£> j)u+iei\(0,c) + e C£l j {J'iJ)ui;{c,oo) 
K 3 

K 3 = (J(0), J) - e l (-+^) c ( J( C ), J) + ( j; jX +lei;(0iC) + 
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+e-^ l{U, + t£l) [e iuc (J(c), J) + < j; J)*.^]. 

At the small values of value e± this expression close to (f5T|) . From it is 
visible that lim e2 _o£2(r) = e~ £lC . 

Unfortunately, to compare the got results to the experiment difficultly, be- 
cause parameters c,a,k,e\ are unknown. But, apparently, through the choice 
of type of distribution p q (u) and selection of his parameters it is possible to 
obtain very good accordance with experimental results. 



6 Comparison with the theory of transport pro- 
cesses by McLennan 

In works [131 EI] the statistical theory of transport processes, based on in- 
troduction of external forces of unpotential behavior, describing influence of 
surroundings on this system, is built. In the Appendix II to work pQ this 
theory is compared with the method of NSO. 

The function of distribution of the complete system f u is described by Li- 
ouville equation 

^ + {f U)Hu } = 0, (32) 

where {, } are classic Poisson brackets, H u = H + H s + U , H is the Hamiltonian 
of the concerned system, H s is Hamiltonian of surroundings, U is Hamiltonian 
of interaction of the system with surroundings. The functions of distributions 
of the concerned system / and surroundings g are accordingly equal 

/ = J fudT s , g = J f u dT, 

where dT s and dT are elements of phase volume of surroundings and concerned 
system. Integrating (j32)) on phase space of surroundings dT s , we get equation 
of motion for /, containing a source. At introduction of function X describing 
correlation of the system with surroundings 

/« = fgx, (33) 

equation for / is written down in a form 

| + {/ ,„ ) + ^ = 0, ,34, 

where 

F a = - [ gX^dr s , (35) 
J 9q a 

q a and p a are coordinates and impulses of the system, in (|34(1 summation up 
is assumed on a. A value F a makes sense "force" presenting the action of 
surroundings on the system. In [431 144j from the physical considering justice 
of expression is assumed 

= ~ J 3s(x)ds, (36) 
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where j s {x is density of entropy flow (including the work accomplished above 
the system) , ds is element of surface limiting the system. For the function of 
distribution in [T] expression coincident with expression for NSO at e — > is 
got. If to accept, that correlation between the system and surroundings results 
in replacement of expression (|36[) on expression 



and Y(x, t) = — J p q (t)dt, after integration on time, conducted at the decision 
of equation for the function of distribution, we will get for the function of 
distribution of the system expressions coincident with expressions of the this 
work. 

Indeed, the function of distribution of time of past life t — t$ of the system 
Pq(u) depends and from properties of the system and from surroundings, as 
X. A function X must depend on time. Then, for example, for expression 
4a) with distribution (| 14[) . (Til))) function — J p q (t)dt related to X(t), to the 
moment of time c decreases, as (-2—) fe , and from a moment c assumes a form 

' ^ a+t > ' 

(■^:) k e~' £C e~ et , e — > after thermodynamic limiting transition and tendency 
(r) — > co. Time c comes forward in a role of time of establishment of stationary 
flows. For an example 4c) with distribution <P5|) a function Y(x,t) decreases 
as e~ eit at t < c and as e c (^-^) e -e2t a t t > c, e 2 ^ at (r) -> oo. Thus, the 
function of distribution p q (u) of time of past life of the system (more precisely, 
integral from her with a reverse sign) can be interpreted and in connection with 
the function of correlation X between the system and surroundings. 

7 The conclusion 

As it is specified in work [15] . existence of time scales and a stream of the 
information from slow degrees of freedom to fast create irreversibility of the 
macroscopical description. The information continuously passes from slow to 
fast degrees of freedom. This stream of the information leads to irreversibil- 
ity. The information thus is not lost, and passes in the form inaccessible to 
research on Markovian level of the description. For example, for the rarefied 
gas the information is transferred from one-partial observables to multipartial 
correlations. In work \§\ values e = l/(r) and p q {u) = eexp{— eu} are ex- 
pressed through the operator of entropy production and, according to results 
|15j , - through a stream of the information from relevant to irrelevant degrees 
of freedom. Introduction in NSO function p q {u) corresponds to specification of 
the description by means of the effective account of communication with irrel- 
evant degrees of freedom. In the present work it is shown, how it is possible 
to spend specification the description of effects of memory within the limits 
of method NSO, more detailed account of influence on evolution of system of 
quickly varying variables through the specified and expanded kind of density 
of function of distribution of time the lived system of a life. 

In many physical problems finiteness of lifetime can be neglected. Then 
e ~ l/(r) — > 0. For example, for a case of evaporation of drops of a liquid 
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it is possible to show [IS], that non-equilibrium characteristics depend from 
exp{y 2 };y = e/(2A2) 1 ' 2 , A2 is the second moment of correlation function of 
the fluxes averaged on quasi-equilibrium distribution. Estimations show, what 
even at the minimum values of lifetime of drops (generally - finite size) and the 
maximum values size y = e/(2X 2 ) 1/2 < 1CT 5 . Therefore finiteness of values (r) 
and e does not influence on behaviour of system and it is possible to consider 
e = 0. However in some situations it is necessary to consider finiteness of 
lifetime (r) and values e > 0. For example, for nanodrops already it is necessary 
to consider effect of finiteness of their lifetime. For lifetime of neutrons in a 
nuclear reactor in work [3] the equation for e = l/(r) which decision leads to 
expression for average lifetime of neutrons which coincides with the so-called 
period of a reactor is received. In work [46] account of finiteness of lifetime of 
neutrons result to correct distribution of neutrons energy. 

Use of distributions ((H) and several more obvious forms of lifetime 
distribution in quality p q {u) leads to a conclusion, that the deviation received by 
means of these distributions values In g(t) from In g zu b(t) is no more e ~ 1/ (r). 
Therefore in expression Q additives to the result received by Zubarev, are 
proportional e. This result corresponds to mathematical results of the theory 
asymptotical phase integration of complex systems [H] according to which 
distribution of lifetime looks like p q {u) — exp{— eu) + Xcpi(u) + X 2 ip2(u) + ... , 
where the parameter of smallness A in our case corresponds to value e ~ l/(r). 
Generally the parameter A can be any. 

For distributions of kind (fl"4)) . (fT9|) . ([23]) . having a various form for different 
times, additives to Zubarev's NSO are distinct from zero and for infinitely large 
systems with infinitely large lifetimes. In the present work it is shown, that 
it is possible, for example, for distributions of lifetime of the system, having a 
various appearance at different stages of evolution of system. Such behaviour 
will be coordinated with known division of process of evolution of system into 
a number of stages [16] . For some distributions the effect of " finite memory" 
when only the limited interval of the past influences on behaviour of system is 
observed. 

Probably, similar results will appear useful, for example, in researches of 
small systems. All greater value is acquired by importance of description of 
the systems in mesoscopical scales. A number of the results following from 
interpretation of NSO and p q (u) as density of lifetime distribution of system [6], 
it is possible to receive from the stochastic theory of storage [33] and theories 
of queues. For example, in [33j the general result that the random variable 
of the period of employment (lifetime) has absolutely continuous distribution 
p q {u) = g(u,x) = xk(u — x,u),u > x > is resulted; g(u,x) = in other 
cases, where k(x,t) is absolutely continuous distribution for value X(t) - input 
to system. 

The form of distribution chosen by Zubarev for lifetime represents limiting 
distribution. The choice of lifetime distribution in NSO is connected with the 
account of influence of the past of system, its physical features, for the present 
moment, for example, with the account only age of system, as in Zubarev's NSO 
[U El [6] at e > 0, or with more detailed characteristic of the past evolution of 
system. The received results are essential in cases when it is impossible to 
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neglect effects of memory when memory time there is not little. The analysis 
of time scales as it is noted in |15) , it is necessary to spend in each problem. 

To determination of type of distribution p q {u) it is possible to apply general 
principles, for example, principle of maximum of entropy, as this it was done 
in the work [47) for determination of type of source in kinetic equation. It was 
marked that exists many possibilities of choice of function p q {u). But by virtue 
of a number of the reasons not all functions can be used in NSO. So, in work [14] 
it was marked that NSO can exist only if corresponding retarded theories are 
memory renormalizablc. In general case offered approach can be interpreted 
as working out in detail of history of evolution of the system, clarification of 
different stages of its conduct. 
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